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Abstract. We use an updated version of SimProp , a Monte Carlo simulation scheme for 
the propagation of ultra-high energy cosmic rays, to compute cosmogenic neutrino fluxes 
expected on Earth in various scenarios. These fluxes are compared with the newly detected 
IceCube events at PeV energies and with recent experimental limits at EeV energies of the 
Pierre Auger Observatory. This comparison allows us to draw some interesting conclusions 
about the source models for ultra-high energy cosmic rays. We will show how the available 
experimental observations are almost at the level of constraining such models, mainly in 
terms of the injected chemical composition and cosmological evolution of sources. The results 
presented here will also be important in the evaluation of the discovery capabilities of the 
future planned ultra-high energy cosmic ray and neutrino observatories. 
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1 Introduction 

Ultra-high energy cosmic rays (UHECRs) are the most energetic particles observed in nature, 

with energies up to a few times 10 20 eV. The first detection of such particles dates back to 

1960 and so far, after several decades of observations, the origin of UHECRs still remains an 

open question. 

The most recent and advanced experiments observing UHECRs are the Pierre Auger 

Observatory [1, 2], by far the largest experimental set-up devoted to the study of UHECRs, 

located in the southern hemisphere in Argentina, and the Telescope Array (TA) experiment 

[3, 4], with 1/10 of the Auger statistics at the highest energies, located in the northern 

hemisphere in the United States. The two experiments adopt a hybrid approach that allows 

events simultaneously observed with a surface detector and a fluorescence detector to be 

collected. The method strongly reduces the uncertainties on energy calibration of the large 

majority of the showers only detected by the surface arrays. 

The experimental efforts made in the field of UHECRs produced a number of significant 

results concerning these particles, such as [5, 6, 7]: (i) UHECRs are charged particles, with 

limits on neutral particles around 10 19 eV at the level of few percent for photons and well 
below for neutrinos [8, 9, 10], (ii) the spectra observed on Earth show a slight flattening at 
energies around 5 x 10 1 * * * * * * * * * * * * * * * * 18 eV (called the ankle), with (iii) a steep suppression at the highest 
energies around 10 2 ° eV [11, 12], 

A key piece of information in the physics of UHECRs is the composition of these ex¬ 
tremely energetic particles. In this respect experimental observations do not converge to the 
same conclusions, with different experiments pointing toward different interpretations. The 
measurements performed by TA are consistent with predictions for a pure proton flux in the 
whole energy range observed, starting from the lowest energies 10 18 eV up to the highest (see 
[13] and references therein), whereas those by Auger signal a richer phenomenology with a 
flux dominated by protons and light nuclei at low energies but becoming progressively heavier 
for E > 3 x 10 18 eV [12, 14], albeit using substantially different analysis modalities. On the 
other hand, TA measurements have wider error bars and a recent analysis [15] has shown 
that they are also compatible with Auger data within their uncertainties. 

During their journey from the source to the observer, UHECRs cover cosmological 
distances and, apart from magnetic fields not considered here, interact with the radiation 
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fields of the astrophysical backgrounds, mainly the cosmic microwave background (CMB) and 
the extragalactic background light (EBL). The propagated energy spectrum of nucleons 1 
is affected almost only by the CMB radiation field and the processes that influence the 
propagation are: (i) pair production and (ii) photo-pion production [19]. On the other hand, 
the propagation of heavier nuclei is affected also by the EBL and the interaction processes 
relevant are: (i) pair production and (ii) photo-disintegration [16, 17, 18, 20, 21]. 

In this work we will consider both scenarios, proton-dominated and mixed composition, 
that are adopted in the interpretation of the most recent experimental data. A pure proton 
composition will be used in the analysis of TA data: under this assumption the observations 
can be described in the framework of the dip model [19, 22], which gives an explanation of 
the ankle observed in the spectrum as an effect of the pair production by protons on the 
CMB, while the strong flux suppression at the highest energies is the effect of the photo-pion 
production process, the so-called Greisen-Zatsepin-Kuzmin (GZK) cut-off [23, 24], 

The Auger data, on the other hand, indicate that the overall proton contribution to the 
spectrum is mostly restricted to the lowest energies < 3 x 10 19 eV well below the photo-pion 
production threshold (~ 6 x 10 19 eV) [25, 26]. The ankle feature has to be interpreted as the 
transition between two different classes of sources, the high-energy one with large metallicity 
while the flux suppression at the highest energies is generated by nuclei photo-disintegration 
and/or the maximum acceleration energy [25, 27, 28]. 

As was first realised by Berezinsky and Zatsepin [29], the proton content of UHECRs, 
particularly at the highest energies, is a crucial quantity that affects the fluxes of secondary 
photons and neutrinos produced by UHECR propagation in the intergalactic space [30, 31]. 
In the present paper, we will focus on the connections between UHECR chemical composi¬ 
tion and source cosmological evolution and the production of secondary neutrinos. In the 
framework of different source models, we will give a comprehensive picture of the expected 
UHECR and neutrino fluxes in comparison with the latest experimental data (previous es¬ 
timates of these quantities by different authors include [32, 33, 34, 35, 36, 37, 38, 39, 40]). 
Particularly relevant are the recent observations performed by the IceCube observatory in 
Antarctica, which detected, for the first time, an extraterrestrial neutrino flux at energies 
around 10 15 eV [41, 42], together with the limits on neutrino fluxes fixed by Auger [43]. 

The calculations presented here are based on SimProp 2 [44, 45], a Monte Carlo (MC) 
propagation code designed to study UHECR propagation and upgraded here to determine 
the fluxes of secondary neutrinos, taking into account all possible production channels. The 
production of secondary neutrinos by the propagation of UHECRs is contributed by the 
overall universe [32], therefore the cosmological evolution of both (i) the sources of UHECRs 
and (ii) the astrophysical backgrounds play an important role. While the impact on the 
expected neutrino fluxes of the latter will be one of the main findings of this work, the poor 
knowledge of the former represents a source of uncertainty that should be quantified. The 
cosmological evolution of the EBL is inferred from (few) observations at different redshifts 
implementing ad hoc models [46, 47, 48]. These models show sizeable differences only at high 
redshift (z > 4), not actually relevant in the determination of UHE nuclei fluxes but affecting 
the production of secondary neutrinos [18]. In the present paper we will use the two recipes 
proposed in [46, 47] and [48] showing the impact of the considered EBL evolution model on 
our findings. 

1 Hereafter discussing freely propagating UHE nucleons we will always refer only to protons because the 
decay time of neutrons is much shorter than all other time scales involved [16, 17, 18]. 

2 SimProp is available upon request to SimProp-dev@aquila.infn.it. 
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The results discussed here, addressing the expected fluxes of secondary neutrinos through 
a detailed analysis of some common UHECR models, has a twofold interest: on one hand, as 
stated above, using the latest observations we can already draw interesting conclusions on the 
sources of UHECRs, on the other hand the study presented should be intended as a bench¬ 
mark computation to assess the discovery capabilities of the next generation experiments in 
both fields of UHECR [49] and neutrino observations [50, 51, 52], 

The paper is structured as follows: in Section 2 we illustrate the physics of the propa¬ 
gation of UHECR through intergalactic space, with a particular emphasis on the processes 
of secondary particles production and on the related new features implemented in SimProp ; 
in Section 3 we illustrate different scenarios for UHECR sources and their output in terms of 
secondary neutrinos, focusing on the experimental limits; finally, in Section 4, we summarise 
and discuss our results. 


2 UHECR propagation and secondary neutrinos production 


Ultra-high energy (UHE) protons or nuclei propagating in the intergalactic space interact 
with photons of the CMB and EBL; these interactions are particularly important not only 
as a mechanism of energy losses for the propagating particles but also being the processes 
responsible for the production of secondary particles such as neutrinos, gamma rays and 
electron-positron pairs 3 . The interaction rate associated to such processes can be written in 
a very general form as [16, 53]: 


1 
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1 
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where T is the Lorentz factor of the particle, <j(d) is the total cross-section associated to the 
particle interactions, d is the background photon energy in the particle rest frame, e( nin is 
the lowest value of d above which the interaction is kinematically possible (threshold), and 
n 7 (e)de is the number per unit volume of background photons with energy between e and 
e + de in the laboratory reference frame. The photon energy in the particle rest frame is 
related to that in the laboratory frame by d = Te(l — cos 0 ), where 0 is the angle between 
the particle and photon momenta (0 < d < 2 Te). 

The interaction processes that involve UHECRs with background photons are pair- 
production and photo-pion production in the case of protons, and also photo-disintegration 
in the case of heavier nuclei. Given the distribution of background photons and the energies 
involved, the propagation of protons is substantially affected only by the CMB radiation 
field, while in the case of nuclei, and only for the photo-disintegration process, also the EBL 
field plays an important role [16, 17]. As we will discuss below, the effect of EBL on proton 
propagation has an important role for the production of secondary neutrinos, but it negligibly 
affects the expected proton flux. 

Let us start from the process of photo-pion production. Nucleons (IV), whether free or 
bound in nuclei, with Lorentz factor T > 10 10 interacting with the CMB photons give rise to 
the photo-pion production process: 


N + 7 —> N + 7T° N + 7 —> N + . (2-2) 

3 In the present paper we will not discuss the case of secondary gamma rays and electron-positron pairs, 

leaving this subject to a possible forthcoming publication. 


- 3 - 




At lower energies T < 10 10 , even if with a lower probability, the same processes can occur 
on the EBL field. In the case of UHE protons propagating in the CMB, the photo-pion 
production process involves a sizeable energy loss producing the so-called GZK cut-off [23, 24], 
a sharp suppression of the flux of protons expected on Earth at E ~ 6 x 10 19 eV, which 
corresponds to the threshold for photo-pion production, that in the nucleon rest frame reads 
4in = m n + m^/2mjv ~ 145 MeV. The photo-pion production cross-section has a complex 
behavior with a number of peaks corresponding to different hadronic resonances, the largest 
one being the A resonance placed at e' = ~ 340 MeV [19]. At energies much larger 

than eA the cross-section sits on an approximately constant value [19]. The photo-pion 
production process holds also for nucleons bound within UHE nuclei, being the interacting 
nucleon ejected from the parent nucleus, but this process is subdominant with respect to 
nucleus photo-disintegration except at extremely high energies [18]. 

UHE nuclei propagating through astrophysical backgrounds can be stripped of one or 
more nucleons by the interactions with CMB and EBL photons, giving rise to the process of 
photo-disintegration: 

(A, Z) + 7 —» {A — n, Z — n') + nN (2.3) 

n ( n') being the number of stripped nucleons (protons). In the nucleus rest frame the energy 
involved in such processes is usually much less than the rest mass of the nucleus itself, 
therefore in the laboratory frame all fragments approximately inherit the same Lorentz factor 
of the parent nucleus, i.e. we can neglect nucleus recoil. The cross-section is dominated by 
a smooth peak, the giant dipole resonance, that appears for photon energies close to the 
threshold (8 MeV ~ e' min < e' < 30 MeV) [20]; in this regime nucleons collectively behave 
as a two fluid system: a proton fluid vibrating against a neutron fluid. The giant dipole 
resonance corresponds to the extraction of one nucleon and it is the dominating process in 
UHE nuclei propagation [16, 17, 20, 21]. At larger energies e' > 30 MeV the quasi-deuteron 
process dominates, in which the photon interacts with one or two nucleons inside the nucleus 
with the extraction of two or more nucleons. This regime corresponds to an almost constant 
cross-section and has a small impact on the propagation of UHE nuclei [16, 17, 20, 21], 

The interaction rate associated to the processes of photo-pion production and photo¬ 
disintegration can be written using Eq. (2.1) specifying the cross-section, the background 
photon density and all relevant kinematical thresholds of the process. 

Protons and nuclei with Lorentz factor T > 10 9 can undergo the process p + 7 —>■ 
p + e + + . The mean free path associated to pair production is relatively short compared 

with all others length scales of UHECR propagation, with a very small amount of energy 
lost by the propagating particle in each interaction. Therefore, as was shown elsewhere 
it is possible to neglect the stochasticity of the process assuming a continuous energy loss 
evolution [19]. 

Particles covering cosmological distances feel the effect of the changes in the background 
universe due to cosmology. The expansion of the universe produces an adiabatic energy loss 
mechanism that continuously reduces the energy of the propagating particles. Assuming 
standard cosmology we can write the energy lost (adiabatically) per unit time by UHECRs 
as 

(“fS") i = ^) = W(l+^m + 0 A (2.4) 

where z is the redshift at time t, Ho ~ 70 km/s/Mpc is the Hubble constant, H m ~ 0.3 is 
the matter density, and Ha ~ 0.7 is the dark energy density. 
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Let us now concentrate on the unstable particles (pions, free neutrons and unstable 
nuclei) produced by the propagation of UHECRs through photo-pion production and photo¬ 
disintegration. In most cases the decay length of such particles is much shorter than all other 
relevant length scales, so these particles decay very soon giving rise to secondary neutrinos. 

There are two processes by which neutrinos can be produced in the propagation of 
UHECRs: 

• the decay of charged pions produced by photo-pion production, —> /u, ± + and 

the subsequent muon decay /u, ± -> eA + /y p (;y /t ) + v e (u e ): 

• the beta decay of neutrons and nuclei produced by photo-disintegration: n p+e~-\-P e , 

{A, Z) ( A , Z — 1) + e + + z^ e , or (A, Z) —> [A, Z + 1) + e + v e . 

These processes produce neutrinos in different energy ranges: in the former the energy of 
each neutrino is around a few percent of that of the parent nucleon, whereas in the latter it 
is less than one part per thousand (in the case of neutron decay, larger for certain unstable 
nuclei). This means that in the interactions with CMB photons, which have a threshold 
around T > 10 10 , neutrinos are produced with energies of the order of 10 18 eV and 10 16 eV 
respectively. Interactions with EBL photons contribute with a much lower probability than 
CMB photons, affecting a small fraction of the propagating protons and nuclei. Neutrinos 
produced through interactions with EBL, characterised by lower thresholds, have energies of 
the order of 10 15 eV in the case of photo-pion production and 10 14 eV in the case of neutron 
decay. 

The results presented in this work are obtained using SimProp v2r2 [45], an extension 
for the study of secondary neutrinos of the SimProp [44] computation scheme. The most 
relevant technical details and a more detailed discussion of the algorithm and numerical code 
implemented can be found in Ref. [45]. The following new features have been implemented 
in SimProp: 

• photo-pion production is taken into account as a stochastic process, for protons as well 
as nuclei, using both the EBL and CMB radiation fields; 

• all allowed isobars for each A, stable or not, can be produced, though unstable ones 
are assumed to decay immediately so only stable isobars for each A are propagated; 

• the decay products of pions, neutrons and unstable nuclei are tracked; 

• several technical improvements were implemented, resulting in a faster execution time; 

A SimProp run consists of N events (for a user-supplied value of N), each consisting 
in the generation of a primary particle and its propagation from the source to the observer 
along with that of any secondary particles produced. All particles are assumed to travel 
rectilinearly (without taking into account the possible presence of intergalactic magnetic 
fields) at the speed of light, so only one coordinate, the redshift z, is used to keep track of 
particle positions. 

In order to test our algorithm, particularly the new module of the code developed to 
compute the flux of secondary neutrinos, we have compared our results with those of a simu¬ 
lation performed by Engel, Seckel and Stanev [33]. In figure 1 we plot the fluxes of secondary 
neutrinos and anti-neutrinos u e , z/ e , z/ /x , z/ /i; where results of [33] are plotted as smooth lines 
while our results through histograms. The two simulations show results that are in quite 
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Figure 1. Fluxes of neutrinos at Earth expected in the EngeUSeckel-Stanev [33] scenario. The 
histograms were obtained via SimProp and the smooth curves refer to the data from Ref. [33]. [Left 
panel] v e flux (blue solid) and v e flux (red dashed). [Right panel] flux (blue solid) and flux (red 
dashed). 


good agreement. The only substantial difference is that our electron antineutrino spectrum 
lacks the second peak at around 10 18 eV. This follows from our choice of neglecting the sub¬ 
dominant interaction channels in which several pions and/or heavier mesons are produced. 
This simplification is observationally irrelevant, because the phenomenon of neutrino oscilla¬ 
tions (not shown here) will mix the flavours so that at Earth all neutrino flavours will show 
equal fluxes and the total antineutrino production at around 10 18 eV is dominated by h /; . 
Neutrino fluxes presented in the forthcoming discussion are always the sum of neutrino and 
anti-neutrino fluxes over all flavours. 

3 UHECR scenarios and cosmogenic neutrino fluxes 

Theoretical models aiming at the explanation of UHECR observations can be distinguished in 
two different scenarios depending on the chemical composition: on one hand the dip model, 
which assumes a pure proton composition, on the other hand models based on a mixed 
composition with the low energy part of the spectrum dominated by protons and the highest 
energy part dominated by nuclei. As we will show below, the flux of secondary neutrinos, 
being extremely sensitive to the chemical composition, can be a powerful tool to solve the 
puzzle. 

Neutrinos produced by the interaction of UHECRs, because of their extremely low 
interaction rate, arrive on Earth unmodified, except for redshift energy losses and flavor 
oscillations, with the overall universe contributing to their flux. This is an important point 
that makes neutrinos a viable probe not only of the chemical composition of UHECRs but also 
of the cosmological evolution of sources that, as we will show below, can be also constrained by 
the neutrino flux observed on Earth. We will consider the case of sources with no cosmological 
evolution, with the same cosmological evolution as that active galactic nuclei (AGN), an 
astrophysical object supposed to play a role in particle acceleration to the highest energies 
[54], and with the cosmological evolution of the star formation rate (SFR). 

The production of UHECRs at the source will be modelled as a power law in the Lorentz 
factor of the injected particles T. We will also account for the maximum acceleration energy 
that sources can provide, T max . This last quantity will be taken to be rigidity dependent, 
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namely given the maximum acceleration energy for protons Emax the same quantity for any 
nuclei species will be E A ax = ZEmax, Z being the atomic number (charge) of the nucleus. 
The spectrum of cosmic rays effectively injected in the intergalactic medium by a collection 
of sources can be affected by the convolution of the spectrum of individual sources and their 
luminosity and/or maximum energy achieved. For instance, under realistic assumptions [55], 
the effective spectrum can have a broken power-law dependence. For soft injection (spectral 
index j g > 2 ), which is mainly of interest in the case of sources injecting only protons and 
light nuclei, we will assume as in [25] an effective injection spectrum of the type 


Qmj(r, z) oc S(z)e r/F; 


r 2 
1 0 


i/r 2 

_r \ “ 79 

W 


r < r 0 
r > r 0 . 


(3.1) 


The value of To depends on the specific adopted model; in the case > 2 for phenomeno¬ 
logical purposes typical values are To — 10 6 -r 10 8 , as this makes the requirements in terms of 
the source luminosity less severe [19]. In this work we use To = 10 8 . Following the approach 
of [25], in the case of sources also injecting heavier nuclei, we are especially interested in hard 
injection, where the slope of the generation spectrum is < 2 ; in this work we do not use 
any spectral break for such sources. 

The term S(z ) in equation (3.1) accounts for the cosmological evolution of sources. 
In the forthcoming discussion we will consider separately three different cases: the case of 
no evolution with S(z ) = 1, the case of the cosmological evolution associated to the SFR 
[38, 56, 57]: 

f (1 + z) 3A z < 1 

Ssfr{z) = { 2 3 - 7 (l + z )-°- 3 1 < ^ < 4 (3.2) 

[ 2 3 - 7 5 3 - 2 (l + z )" 3 ' 5 z>4 

and the case of the evolution typical of the AGN [38, 58, 59] 

({l+zf'° z <1.7 

Sagk(z) = (1 + 1.7 ) 5,0 1.7 <z< 2.7 (3.3) 

[ (1 + 1.7) 5 -°10 (2 - 7 "^ z>2.7 


All computations presented here are performed under the assumption of a homogenous 
distribution of sources. This assumption does not affect the expected neutrino spectra be¬ 
cause in the case of neutrinos the overall universe, up to the maximum redshift, contributes 
to the fluxes and possible flux variations due to a local inhomogeneity in source distribu¬ 
tion gives a negligible contribution to the total flux. We also fix a maximum redshift of the 
sources z max = 10, which is the typical redshift of the first stars (pop III) [60]. In any case 
the expected fluxes of primary and secondary particles are almost independent of z max if 
Zmax 3. 

Once produced at cosmological distances neutrinos travel toward the observer almost 
freely, the opacity of the universe to neutrinos being relevant only at the redshifts z > ID 
[61, 62]. Therefore, given the assumptions discussed above, in our computations we have 
neglected any effect due to neutrino propagation apart from the adiabatic energy losses due 
to the expansion of the universe as discussed in the previous section. 

In the following subsections we will consider the two cases of UHECRs consisting of a 
pure proton composition or a mixed composition, assessing their consequences in terms of 
secondary neutrino production. 
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Figure 2. [Left Panel] Fluxes of protons expected at Earth in proton-only scenarios with various 
models for the cosmological evolution of sources (solid red: no evolution; dashed green: SFR evolution; 
dot-dashed blue: AGN evolution [38]) normalized to TA data (purple squares); Auger data (olive 
disks) and KASCADE-Grande data [63] (blue triangles) are also shown for comparison. [Right Panel] 
Fluxes of neutrinos in the same scenarios (same color code; from bottom to top: no evolution, SFR, 
and AGN), with colored bands showing the difference between the Stecker (solid) [46, 47] and Kneiske 
(dashed) [48] EBL models; thin solid lines are neutrino fluxes obtained taking into account only the 
interaction with CMB photons. 

3.1 The dip model 

The dip model is based on the assumption that UHECRs are composed of only protons, 
therefore in this case we will restrict our analysis to the TA data. In order to describe the 
overall observed spectrum in terms of a single proton component we need to assume a rather 
steep injection spectrum with > 2.3 [19, 22], The actual value of the injection spectrum 
depends on the assumptions about the cosmological evolution of sources: increasing the 
co-moving density of sources at large redshift the required value of decreases [19, 22]. 

In figure 2 we consider the three different scenarios for cosmological evolution of sources 
mentioned above, namely: no evolution, SFR and AGN evolution, assuming that contributing 
sources are restricted to redshifts z < z max = 10. In order to reproduce the TA observations 
shown in the left panel in figure 2, the injection power law index r y g needs to be changed 
according to the cosmological evolution: we used = 2.6, 2.5, and 2.4 respectively in 
the three cases of cosmological evolution considered. The maximum acceleration energy 
used to compute fluxes in figure 2 is Emax = 10 22 eV. The source emissivities, i.e. the 
energy of UHECRs emitted per unit volume and time, required to reproduce TA data are 
£o ~ 1-5 x 10 46 , 6 X 10 45 and 3.5 x 10 45 erg/Mpc 3 /yr at z = 0 respectively in the case of 
7 g = 2.6, 2.5, and 2.4, assuming To = 10 8 (see equation (3.1)). 

Results presented in figure 2 are quite interesting, showing that models with cosmolog¬ 
ical evolution of sources stronger than the case of AGN predict a total neutrino flux (right 
panel of figure 2) in excess of the experimental limits fixed by Auger and the observations 
of IceCube. Moreover, as predicted in [54], the AGN hypothesis seems to reproduce quite 
well TA observations in terms of the UHE flux suppression, where our calculations result in 
a neutrino flux at a level which would have been detected by both IceCube and Auger. 
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3.2 Mixed composition 

Let us consider here the case of a mixed composition, which is more suitable to describe the 
observations of Auger; in this analysis we will refer only to Auger data on both spectrum 
and chemical composition. 

As shown by different authors [25, 26, 64] scenarios with mixed composition have pe¬ 
culiar characteristics that, in order to describe Auger observations, imply different source 
families (see [65, 66 ] for models with a single type of source). At the highest energies 
(E > 5 X 10 18 eV) Auger observes predominantly heavy elements, with A > 12, while at 
the lowest energies, as a matter of fact, all UHECR experiments observe a light composition 
dominated by protons. To reproduce such observations two classes of sources are needed 
[25, 26, 64]: one that injects only light elements with low maximum energy and steep in¬ 
jection and one that injects also heavy elements with higher maximum energy and a flatter 
injection. 

In this context, while the extragalactic nature of sources that contribute also heavy 
elements cannot be questioned, sources providing only light elements can be both galactic or 
extragalactic, though the galactic hypothesis poses severe problems with both the limits on 
observed anisotropy and the (galactic) acceleration mechanism [25] . In this case the produc¬ 
tion of secondary neutrinos will be connected only to the propagation of UHECRs injected 
by extragalactic sources and will result in a particularly low flux, far below the detection 
threshold of any experiment. For a different approach with galactic sources providing the 
low energy light component and secondary EeV neutrinos produced directly at the source 
the reader can refer to [67, 68 ]. 

In the present paper, in order to study the production of cosmogenic neutrinos in the 
framework of the Auger observations, we will consider the case discussed in [25], but with 
slightly retuned spectral parameters in order to account for the different approximations used 
in SimProp and those used in [25] and the differences between the 2009 Auger data used in [25] 
and the 2013 data used here. The model is composed of two different classes of extragalactic 
cosmic ray sources distributed uniformly in a comoving volume. The first class injects (75%) 
protons and (25%) helium nuclei with = 2.6, 2.5, and 2.4, depending on the cosmological 
evolution as in the previous section, with a maximum energy E max = Z x 2 x 10 18 eV and a 
source emissivity Cq ~ 1.5 x 10 46 , 5.5 x 10 45 and 2 x 10 45 erg/Mpc 3 /yr at z = 0 respectively 
in the three cases = 2.6, 2.5 and 2.4, assuming To = 10 8 (see equation (3.1)). The second 
class also injects heavier nuclei, with 35% protons, 30% helium, 25% CNO and 10% MgAlSi 
(these fractions being defined at fixed energy per nucleus equal to 10 18 eV), with 7 = 1.0, 
E m ax = Z x 6 x 10 18 eV and a source emissivity Cq ~ 5.5 x 10 44 erg/Mpc 3 /yr (above 
10 7 GeV/n). 

It is interesting to note here that in the case of nuclei the cosmological evolution of 
sources has a larger impact on the expected fluxes than in the case of a pure proton com¬ 
position. This can be understood taking into account that cosmological evolution implies an 
increased source density at high redshift, through terms of the type (1 + z) m (see above), 
that corresponds to an increased effect of the photo-disintegration interaction. This fact 
produces an increased number of secondary light cosmic rays, particularly at intermediate 
energies E ~ 10 18 ' 5 eV, spoiling the agreement with both flux and chemical composition 
observed by Auger. Therefore sources providing heavy elements cannot show any evolution 
with redshift. This result is actually quite interesting, showing the capability of the Auger 
data on flux and chemical composition in constraining also the cosmological characteristics 
of the possible UHECR sources. 
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Figure 3. Fluxes of UHECRs (left) and neutrinos (right) expected at Earth in the two-class model. 
UHECR fluxes are grouped according to their mass number at Earth as labeled. Olive disks in the 
left panel are Auger data. TA data (purple squares) and KASCADE-Grande data [63] (blue triangles) 
are also shown for comparison. Neutrino fluxes are shown (from bottom to top) for the homogeneous 
(red), SFR (green) and AGN (blue) evolution of the steep (light) component, assuming the Stecker 
(solid) and Kneiske (dashed) EBL evolution models. Fluxes of protons and He are summed over the 
two different classes of sources. For the injection chemical composition and source emissivities see 
text. 




Figure 4. Mean (left) and standard deviation (right) of X rnax predicted by EPOS-LHC (solid), 
QGSJet 11-04 (dashed) and SIBYLL 2.1 (dot-dash) for pure protons (top, blue), pure iron (bottom, 
red), and our two-class model (in the no-evolution scenario), along with Auger data (lines: statistical 
errors; brackets: systematic errors). 

Using two classes of sources, with the parameters discussed above, and assuming no 
cosmological evolution of those providing also heavy elements, we obtain the fluxes at Earth of 
UHECRs and neutrinos shown in figure 3, left and right panel respectively. The corresponding 
chemical composition is shown in figure 4. 

The results presented in left panel of figure 3 and in figure 4 confirm the capability of 
models based on two classes of sources of reproducing the Auger observations, as previously 
found in [25, 26, 64], 

The corresponding flux of secondary neutrinos is extremely low at the highest energies 
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as follows from the right panel of figure 3. This is a direct consequence of the chemical 
composition that at the highest energies is dominated by heavy nuclei. In this case the 
production of secondary neutrinos, being mainly connected with the photo-pion production 
process, involves only nucleons inside nuclei with energy E/A. Therefore, assuming a nucleus 
maximum energy at the level of E max = Z x 5 x 10 18 eV, as in figure 3, implies that only 
the tails of the cosmic rays spectra contribute to the production of secondary neutrinos with 
a substantial reduction of their expected flux. Only the light component of UHECRs, being 
proton dominated, will provide a sizeable amount of secondary neutrinos. Therefore, given 
the constraints introduced by the Auger data, we can conclude that only the PeV component 
of cosmogenic neutrinos can be detected, with an expected flux close to the IceCube sensitivity 
assuming a moderate (SFR) cosmological evolution for the light component of UHECRs and 
in excess of it for a strong (AGN) evolution, whereas EeV neutrino fluxes are way below the 
Auger sensitivity. 

4 Conclusions 

Cosmogenic neutrinos mainly originate from the decay of pions (and muons) produced in 
photo-hadronic interactions of protons and heavy nuclei with CMB and EBL backgrounds. 
These processes are efficient only in the case of protons, while in the case of heavier nuclei 
photo-hadronic interactions are significantly suppressed. Therefore the bulk of secondary 
neutrinos is produced by the lighter components of UHECRs. 

There is a solid consensus about the light composition of UHECRs in the low energy 
part of the observed spectrum. This assures a flux of cosmogenic neutrinos in the PeV energy 
region, produced by the proton’s photo-pion production process on the EBL photons. Our 
remaining ignorance in predicting this neutrino flux is due to the cosmological evolution of 
cosmic ray sources and, to a lesser extent, to the uncertainties in the EBL evolution. Models 
with a strong cosmological evolution, as in the case of AGN, produce a flux of cosmogenic 
neutrinos almost at the level of the IceCube observations in the PeV region. This fact enables 
a partial constraining of cosmological evolution of sources, disfavouring models with a too 
strong evolution: S(z) oc (1 + z) m with m > 3.5. 

At the highest (EeV) energies the flux of cosmogenic neutrinos mainly originates from 
the highest energy tail of the UHECR spectrum. If the most energetic cosmic rays, at 
around 10 20 eV, are mainly protons, then a sizeable flux of cosmogenic EeV neutrinos will 
be produced. In this case, as before, the experimental limits on EeV neutrinos can be used 
to constrain the cosmological evolution of sources. For instance, Auger limits on neutrino 
fluxes already disfavour too strong evolution models, i.e. models with m > 5, and moderate 
cosmological evolution, i.e. models with m < 5, can be tested by future UHE neutrino 
detectors. Conversely, if the highest energy tail of UHECRs is composed mainly of heavy 
nuclei, as in models reproducing Auger data on spectrum and chemical composition, then 
the flux of cosmogenic EeV neutrinos is far below the detection threshold of any running or 
planned detector. 
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